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I. INTRODUCTION 

! This paper contains a study of some technical aspects associated with multiparty Bell inequalities. Multiparty Bell 

■ inequalities are recently gaining in importance, because security, or performance, of many of quantum communication 
' schemes (like multiparty key (secret) sharing [l|, quantum communication complexity problems can be measured 

^ . with Bell inequalities of some form f3|. Bell inequalities show the limit of achievable correlations between many 
O ' parties, if they share only classical means of communication and/or computation. If violated by certain predictions 
for a certain entangled state, they show that this state, without any further manipulations is a good resource for 
some quantum informational tasks. To put is short, in the standard context Bell inequalities show that if there are 
N separated parties, each equipped with a classical supercomputer, who share common data, and perhaps programs, 
^-H ' but otherwise are not allowed to communicate, the locally run programs cannot simulate some phenomena associated 
I> ' with spatially separated N local measurements, with locally decided settings, on N entangled particles, when each 

■ party is asked to predict (compute) result of a measurement for a different particle. 

op ' Bell's theorem was formulated for two particles It took a quarter of century to realize that for more than 
two particles the situation is much more interesting. The discovery of Greenberger-Horne-Zeilinger correlations Q 
^ immediately led to various generalizations of the multipartite problem. This gave birth to the quest of finding Bell 
\^ , inequalities for the new problems. Mermin was first to produce series of inequalities for arbitrary many particles, 

■ involving dichotomic observables, and allowing each observer to choose between two settings [Gj. A complementary 
series of inequalities was introduced by Ardehali ^ . In the next step Belinskii and Klyshko gave series of two settings 
inequalities, which contained the tight inequalities of Mermin and Ardehali [§]. Finally the full set of tight two 
setting per observer, N party Bell inequalities for dichotomic observables was found independently in [9|, [l^, UH- 

■ All these series of inequalities are a generalization of the CHSH ones (T^]. Such inequalities involve only N party 
d correlation functions. The aim of this paper is to go one step further, and to produce series of tight CHSH-type Bell 
^ inequalities for N observers, dichotomic observables, and three possible setting for each observer. Recently series of 

tight inequalities, which do not form a complete set, were found for the case when k-th partner can choose between 
^ 2*^^^ settings (see [l^ [13, Ell)- We shall not follow this approach here, as it seems to be incapable of generating the 
• >— j full set of Bell inequalities. Also the methods of ^16*1 and T?] will be not discussed here, as they do not follow from 
the analysis of the geometry of the polytope of local realistic models (for this concept see [l^l ) . 
^ I Only the proof for the necessary and sufficient condition for the two-observer (3 x 3) case will be shown in full detail. 
- - ' As the reader will see, the generalization to more parties, that is e.g. to 3 x 3 • • • x 3 problems, is straightforward. In 
all these cases we shall be concerned with the generalization of the inequalities of the CHSH type. Such inequalities 
involve only N party correlation functions. However, a remark will be made on how to extend the validity of the 
presented results to CH-type inequalities [l^ (which involve also correlations of lower rank). We shall study only 
experiments involving dichotomic observables. 

We have chosen the three settings case, with two valued observables, because it is the simplest example for which 
the method used here can lead to new results. Also the specific trait of three settings inequalities is that the observers, 
in the case of qubits, can use a full set of mutually unbiased observables for their measurements. The aim of this paper 
is basically a presentation of full details of the derivation of the form of the coefficients in three setting inequalities. 
Explicit forms of inequalities that can be obtained using this formalism will be presented elsewhere. 

Finally note that the results of Garg [2^ are equivalent to the present ones, for the two observer case. However, 
the method presented here is easily generalizable to the case of three or more particles, and therefore extends to an 
unknown territory. Some results specific results concerning three-particle inequalities are presented in [21||. 
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II. DEFINITION OF 3 X 3 PROBLEM 

We shall be interested only in the conditions of existence of local realistic models of the two particle correlation 
functions, forgetting totally the lower order (local) processes. A proper symmctrization of a solution of the present 
problem, leads always to a solution of a problem in which one assumes that lower order (one particle) averages vanish. 
We allow each observer to choose between three settings. We will search for the positivity condition for the unphysical 
(hidden) probability distribution involving results for all settings: 

P{m^^\m^^\mf\mf\m^^\m^^^ (1) 

where m-"'-* is the local realistic prediction for measurement result if the observer i chooses to measure a dichotomic 
observable number m. This distribution must properly reproduce all values of two particle correlation function for 
the full set of 3 x 3 experiments. 

A. Structure of the polytope 

The classical local realistic model of the two particle correlations is essentially a probabilistic mixture of deterministic 
"orders" that are sent by the source to the detection regions. The orders sent to one of the observation stations, in 
the case of three possible settings, must contain the information which value of the local dichotomic observable is to 
result during the measurement. Therefore such an order, for the first observer, can be put in the following way: 

(m;,mj,m2) = (±1,±1,±1). (2) 

We have 2^ distinct orders, that can be sent to one of the observation stations. 

If one is interested only in local realistic models of the two-particle (highest order) correlation functions, then the 
optimal way to describe the deterministic orders is by tensor products of the orders sent to each of the observers. 
Thus the deterministic orders are now: 

{m^,ml,ml) ^ {rr4,ml,ml) = (±1, ±1, ±1) O (±1, ±1, ±1). (3) 

Note however that, the description has an in-built simplification, because one cannot distinguish between 

(m?,mi,mi) (g) (-TO2,-TO2)-"i2) (4) 

and 

{-ml,-ml,-ml) 1^ {m^,ml,ml). (5) 

A local and realistic model of the correlation functions i^mnj) ^i) withn^ = 0, 1, 2, exist if the "tensor" built out of 
these values of the correlation function, E, is reproduced by the model of the following structure: 

PC^(") ^(1) ^(2) (0) (1) (2)^ 

X {m\, m}, m\) ® (m2, rri^)- 

(6) 

As we shall see, the tensor product structure of the full set of correlation functions is a very important property, 
which can be used to simplify the derivations. This one of the basic concepts, around which this paper is formulated. 

III. PEDAGOGICAL INTERLUDE: STRUCTURE OF POLYTOPE OF MODELS FOR JUST ONE SIDE 

In this subsection are presented basic the traits of realistic models for results of experiments on just one side (say 
for Alice). The full set of such models, as we shall see, forms a three dimensional cube. The aim of this section is to 
get some intuition, before one moves to the much more complicated case of two separated observers, and the problem 
of the local realistic models for their observations. An impatient reader can move right away to the next section. 

If one is interested only in the events on one side, say Alice's, of the experiment, then all possible realistic models of 
the average values obtained in the three possible experiments (settings) are convex (i.e., probabilistic) combinations 
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of the local orders Q. Therefore any realistic model lays within the interior of a cube, with eight vertices defined by 
®. 

Now let us understand some of the traits of this cube. First notice that the six faces contain the following vertices: 
the first face 

(l,a,5), (7) 

with a, 6 = ±1, the second face, on the opposite side 

(-l,a,6), (8) 

the third one 

(a, 1,6), (9) 

the fourth 

(a, -1,6), (10) 

fifth 

(a, 6,1), (11) 

and the last one 

(a, 6,-1). (12) 

Always three arbitrary vertices of one face form a complete set in (a non-orthonormal basis). The missing fourth 
vertex from the given face, once one chooses the three vertex basis, is a linear combination of the basis vertices, with 
the expansion coefficients of modulus of one, and the sum of the coefficients equal to one too. For example take the 
first face. The four vertices satisfy the following simple relation: 



a6(l,a,6) = 0. (13) 



a,b=±l 



But one does not have to manipulate this relation in order to make sure that the italicized thesis is correct. Simply 
there is a well known fact from analytic geometry that makes this relation obvious: A hyperplane in an N dimensional 
real space is defined by N linearly independent points that belong to it, and all other points belonging to this hyperplane 
are linear combinations of these N points, with coefficients that have the property that their sum is equal to one. That 
is if fi, with i = 1, N, belong to a hyperplane and are linearly independent, then any other point of the hyperplane 
r is expressible by r = tifi, with ti — I. 

Just to build some intuitive background for the more complicated case to be considered in the sequel, please notice 
the following simple facts about the cube. 

If we take at random three vertices, then these may belong to one face. If this is so, there is a fourth one which 
belongs to this face. Consequently, the plane defined by the three vertices contains a face of the (convex) polytope. 
There are six such planes, and therefore six inequalities (defining half spaces), that have the property that the points 
satisfying all these inequalities form our cube (these could be called local Bell inequalities). 

If the three randomly chosen vertices do not belong to a single face, then they define a plane that cuts through the 
cube. This in turn means that there is a pair of vertices v and —v, which has the property that the vertices are in 
the interiors of the two half spaces defined by the plane, each one in a different one. 

Next, a trivial but important observation. A random triplet of vertices does not have to form a complete set. In 
such a case, there is a vertex v inexpressible by a linear combination of the trio, and the vertex —v has the same 
property. 

Finally, it is easy to see that the relation (fT5|) has the following corollary. If one builds a (nonorthonormalized) 
basis out of three vertices of the cube, then all vertices of the cube are linear combinations of the basis vertices with 
coefficients equal to ±1, or zero. As a matter of fact either three of the coefficients have modulus equal to one, or just 
one of them. 

As we shall see below, the polytope of local realistic models for two-observer experiments has many traits which are 
generalizations of the ones presented above for the cube. These traits can be used to determine the full set of (tight) 
Bell inequalities satisfied by the local realistic models. 
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IV. BACK TO THE CORRELATION POLYTOPE 



In this section the general properties of the faces of the correlation polytope will be shown. These properties define 
the hypcrplancs which contain the faces of the polytope. Therefore they define the Bell inequalities (the full set of 
tight Bell inequalities). 

Now we work in a 9 dimensional real space of vectors. The correlation polytope has 32 vertices. Each vertex v 

has its twin —v. Thus the set of vertices can always be split into two disjoint subsets of 16, such that if v belongs 
to one of them —v belongs to the other one. Each such a subset contains subsets of nine vertices which form a 
non-orthonormalized basis. The question is what is the relation of the hyperplane, defined by the nine vertices, with 
the polytope. Does it contain a face of it, or does it cut through it? In the first case the equation of the hyperplane, 
in an obvious way, leads to a Bell inequality. 

Fact 1. Let us take a random set of nine vertices, U2,...,W9, which forms a complete non-orthonormal basis. 
Assume that this set has a property that 7 more vertices are linear combinations of the nine basis ones, with the sum 
of all expansion coefficients equal to +1. Thus, because of the property of the hyperplane, all 16 vertices are in the 
same hyperplane. Thesis: no other vertex is in the hyperplane, and all those other vertices are on one side of the 
hyperplane (i.e., in the interior of a single half space, out of the two which are defined by the splitting of the space 
by the hyperplane). 

Before we start with the proof, let us have a look at a well known fact from analytic geometry: The equation of 
a hyperplane in a N dimensional space, defined by N linearly independent points w;^, ... , is given by following 
relation. The point x belongs to the hyperplane if 



Z)[a;,w^,t«^, ...^w^] = det 



1 Xi X2 



1 JV JV 
1 W2 



Xn 



w 



N 



w 



JV 



JV -I 



= 0. 



(14) 



D[x, ,w'^, w^] stands for the determinant of a (iV + 1) x (A^ + 1) matrix with the first row beginning with 1, and 
then containing all N coordinates of the point x, that is xi,..., xn, the fc-th row, k = 2,3, ...N, also beginning with 
1, and then containing the coordinates of the point w'' (placed in the same order as those of x). Note that the first 
column of this matrix has all entries equal to 1. The two half-spaces, with the hyperplane as their joint subset, are 



defined by D[x, ,w'^ , w'^] < and D[x, , w'^, ... 
Let us now move to the proof of Fact 1. The vertex, 



„JV 



]>0. 

not belonging to the hyperplane, which contains 16 vertces, 



may be either equal to some —Vk, or is a linear combination X^j^^ Uvi, with ^^^i U = —1. In the first case, one can, 
assume first that = —v^. Then D\—V\,V\,v-2,...,v^ can be expanded using the minors associated with the first 
column. Therefore it reads 



d\px,Vi, Vg] - d[-Vi,V2, Wg] + d[~Vi,Vi,V3, V4, Vg] 
.... - d[-Vi,Vi,V2,—,Vs], 



(15) 



where d[vi,V2,V3, ...,vg] is the determinant of a 9 x 9 matrix with rows defined by the vectors vi, ... vg. Of course, 
such a determinant does not vanish only if all rows of the matrix are linearly independent. Thus 



Similarly, if = —v^, then 



Generally, 



D[-Vi,Vi,V2, ■..,Vg] = 2d[vi,V2, ...,Vg]. 
D[-V2,Vi,V2, ...,Vg] = 2d[Vi,V2, ....,Vg]. 
D[-Vk, Vi,V2, ...,Vg] = 2d[vi,V2, ....,Vg]. 



(16) 
(17) 
(18) 



That is, all points —Vk are in the interior of only one of the two half-spaces. Let us therefore move to the second case 

of = x;i=i tiVi, with x;^! = 



D[Yfi^ltiVi,Vi,V2,...,Vg] = d[Vi,V2,::,V9] - tid[vi,V2, ...,Vg] 
+t2d[v2,Vi,Vs,V4, ...,Vg] + ... - tgd[vg , Vi , V2 , Vg]- 



(19) 
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But, since the opposite vertex —v^ = — X]i=i ^i'^i: by our current assumption, belongs to the hyperplane, one has 

= i'i-X^Ll ti'^i^Vi,V2, ..^Vg] = d[vi,V2, Wg] - ( - tid[vi,V2, ■■■,Vg] 
+t2d[v2,Vi,V3,V4,, ...jWg] + ... - tgd[vg,Vi,V2, ■■■^Vs])- 

(20) 

therefore again we have 

D[vj;,vi,V2, ...jWg] = 2d[vi,V2, ....^vq]. (21) 

QED. 

Fact 2. If vertices Vr and —Vr belong to a hyperplane, then the hyperplane contains the central point of the 
polytope, and thus it cuts through the polytope. 

Before we proceed further, it is good to notice the following simple property of the vertices of the polytope. 

Fact 3. Every vertex Vq is a linear combination of nine vertices Vk, with fc = 1,2, ...,9, which form a complete 
non-orthonormal basis (this is obvious), Vq — Yl\=i tk^k- Take any vertex outside of this set of nine, i.e. certain Vx, 
such that Vx ^ iw^. Then, the expansion cocfRcients of Vx in terms of the basis vertices have their moduli either 
equal to zero, or to a natural number. |,22] Note that this is an extension of one of the properties of the cube. 

Proof. Take nine vertices of the form (l,a',6') (g) {l,c',d') with a',b',c',d' = ±1 , but with the following cases 
excluded a' = b' = 1 and c' = d' ~ 1. The prime is to indicate that not all values ±1 are allowed. Since (1, 1, 1) — 
~S(a' f)')#(i 1) '^'^'(-'^' ^')' the immediate consequence is that (1,1,1) x (l,c',d') have the required property. So 
do (1, a', 6') (g) (1, 1, 1). Finally (1, 1, 1) x (1, 1, 1) = E(a' Mmi.i),{c' a'b'c'd'il, a', b') ® (1, c', d') has also the 
required property. All other vertices are the one listed above with an overall minus sign in front of them, thus 
for them the property holds too. Thus in such a case all expansion coefficients have the property: \tk\ — or 1. 
Take an other set of nine vertices, which form a complete set. In the simplest case, such a new basis can also be 
a "product basis". E.g., it can be defined by, say, excluding cases like, e.g. a" — — 1 and c" — —d" — 1, 
and assigning to the chosen vertices specific signs, i.e. the basis would now be CT(a" , 5" , c" , c?")(l, a" , 6")(1, c" , (i")for 
(a", 6") ^ (1,1) and (c",d") ^ (1,-1), where (T{a" ^b'',^' ,d") ±1 is a "sign" function. This means that with 
respect to the basis, (l,a",6") (l,c",(i"), for which the required property obviously holds, some of the expansion 
coefficients would change their sign, but not the modulus. That is, they again are equal to ±1 or 0. Further, one may 
have a set of nine linearly independent vertices, which cannot be expressed as a set of tensor products of two three 
dimensional bases. E.g., this could be the first basis used in this proof with, say, (1, —1, —1) (g) (1, —1, —1) replaced 
by (1, —1, —1) (g) (1, 1, 1). Therefore now (1, —1, —1) (1, —1, —1) is outside the basis. But it can be expressed as 
(1, —1, —1) (8) [—(1, 1, 1) + (1, —1, 1) + (1, 1, —1)], i.e. the expansion coefficients are of modulus or 1. Other vectors, 
which did not belong to the old basis, and do not belong to the new one, have to be re-expressed in a similar way, 
whenever they were a linear combination which included (1,— 1,— 1)(K)(1,— 1,— 1). It is therefore obvious that their 
new expansion coefficients can only: stay unchanged, or change their values by ±1. Finally, if one considers one more 
replacement in the basis, this again can change the expansion coefficients by only ±1. Therefore, their moduli must 
be equal to 0, 1,2, .... QED. 

Fact 4- Assume now that we choose at random nine vertices, Vk, with k = 1,2, ...,9, and they turn out to form 
a complete non-orthonormal basis. No pair Vr and —v^ belongs to the hyperplane defined by the vertices. This is 
due to the fact, that if the first vertex belongs to the hyperplane, the sum of its expansion coefficients is 1, and thus 
for the opposite vertex it must be —1, i.e. it is out of the hyperplane. Therefore if Vi belongs to the hyperplane, 
then —Vi is out of it. As it was shown earlier, such hyperplane can contain maximally 16 vertices. If it contains 15 
or less vertices, there must be a pair Vx and —Vx which does not belong to it. We shall show that in such a case 
the hyperplane cuts the interval defined by —Vx and Vx into two non-trivial intervals, and therefore it cannot be a 
hyperplane that contains a face of the polytope. That is equivalent to the statement, that the hyperplane cuts through 
the polytope, and inequalities associated with it cannot be used to define the interior of the polytope (i.e., they are 
not Bell inequalities). 

Proof. Let us calculate the value of D[± ti^ii ^^ij ''^2, vg]. If one expands this determinant into determinants 
of minors associated with the first column, one gets 

-D[±Z]i=l UVi,Vi,V2, --^Vg] = d[vi,V2, ■■■■,Vg] ± ( - tid[vi,V2, ■■■,Vg] 
+t2d[v2,Vi,V3,V4,, ...,Wg] + ... - tgd[vg,Vi,V2, -.^Vs]) 

= d[vl,V2,....,Vg]ilTE^t^)■ (22) 

But since it was assumed that neither Vx = Y^^=i ^i^i ^'^^ belongs to the hyperplane defined by u^'s, one must 
have I Y^i=i 7^ 1- It ^ill be shown that | ti\ > 3. 
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First, notice that either Vx or —Vx belongs to the hyperplane Ho containing all vertices of the form (1, a, b) (8) (1, c, d) 
(let us call this set Fg). Next consider a randomly chosen set of any nine linearly independent vertices, B. It can 
be always constructed in the following way. We choose in Fq a set Bq, which contains nine linearly independent 
vertices, and then choose a random subset of these, and we multiply them by —1. I.e., the set B is of the form 
a{a,b, c, d)(l, a,b) ® (l,c, c?), with only a subset. A, of a,d,c,d allowed, and a{a,b, c, d) ~ ±1. If, say, Vx belongs 
to Ho, then the sum of its expansion coefficients, tabcd , with respect to the basis Bo, given by (l,a, 6) ® (l,c, d), 
with a,b,c,d in A, X](a 6 c (i)eA ^"hcrf, is equal to one, and a similar sum for —Vx must give —1. However, the 
coefficients with respect to B, denoted as i^^cs' ^'^^ therefore given by i^fccd ~ o-{0',b,c,d)tabcd- Therefore, by Fact 3, 
t'abcd - tabcd = 0, ±2, ±4, and thus the value of 'E{a,b,c,d)&At'abcd can differ from 'E{aAc,d)eAtabcd by 0, ±2, ±4, .... 
I.e. either J2(a b c d)eAt'abcd ~ which is impossible by our assumption that both zLvx do not belong to H, or 
\^{a b c d)£A t'abcdl — 3- This implies that the value of fF^ . which is non-zero for both vertices Vx and —Vx, is negative 
for one of them and positive for the other. That is, they are in the interiors of the two different halfspaces separated 
by H. I.e., H cannot contain a face of the polytope. It cuts through it. QED 

The conclusion of this section is the following one. All faces of the polytope contain 16 vertices. 

A. Bell inequalities 

A straightforward application of the formula , and of the other above results implies that the tight Bell inequal- 
ities are obtainable by the following procedure 23] . If we find nine vertices, which generate a hyperplane containing 
a given face of the polytope, and arrange them is a sequence in such a way that d[vi, V2, ■■■■,vq] > 0, then the (tight) 
Bell inequality reads 

D[E,vi,V2,...,V9]>0, (23) 

where E — (Eqq, Eqi, Eq2, Eiq, En, E12, E20, E21, £22)- The tensor products defining the vertices, v = {x,y,z) ® 
(w, V, t), are written in a convention consistent with the one for E, that is u = (xw, xv, xt, yw, yv, yt, zw, zv, zt). 

V. BELL INEQUALITIES: THE EXPLICIT FORM 

In order to show how explicitly such inequalities look like, and to have an easy method to pinpoint them, one should 
still: 

• Show how to find which of the sets of 16 vertices belong to one face. 

A. Basis sets of nine 

If we have 16 vertices of one face, then a basis generating a hyperplane containing all these vertices is their 
arbitrary subset of nine linearly independent vertices. Therefore we can always choose a product basis of the form 
a{a, b, c, d){l, a, c) (g) (1, c, d), with (a, b) ^ (1, 1) and (c, d) ^ (1, 1), where a{a, b, c, d) = ±1. The only question left is, 
for which sign functions, (T(a, b, c, d), such a set of nine generates a hyperplane containing a face. 

B. How to find the sets of 16 vertices that generate Bell inequalities (polytope faces) 

Definitely the set (1, a, 6)® (1, c, d), with a, b,c,d = ±1, is in one face of the polytope. Any other set of 16 vertices that 
does not contain oppositely pointing pairs of vertices, has vertices of the following form a{a, b, c, d)(l, a, b) ® (1, c, d), 
with a, b,c,d = ±1, and a specific "sign" ' function cr(a, 6, c, d) = ±1. The problem is to decide whether the given set 
belongs to a single face of the polytope. 

We shall show that if and only if a{a,b,c,d) — S{a,b,c,d) where S{a,b,c,d) is such that J2a b'^^^i'^^^^'^^^) ~ ^ 
and 'Y^^ ^cdS{a,b,c,d) = 0, then the full set of 16 vertices is in one hyperplane (which in turn means that they 
belong to one face). Of course the two conditions warrant also that one has, e.g., ^abcS{a,b, c, d) — and 

bed o-ticdS{a, b, c, d) = 0. Therefore S{a, b, c, d) must have the following form: 



S{a, b,c,d) ^ X + Aa + Bb + Cc + Dd + Eac + Fad + Gbc + Hbd, 



(24) 
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where X,A,....,H are constants. 

Proof: For a basis wc can choose a specific basis set, that is S(a' , b' , c' , d')(l, a', b') ® (1, c', d'), which excludes the 
seven verices S{1, f , c, d){l, f , f ) ® (f , c, d), and S{a, b, 1, f )(f , a, b) ® (1, f , 1). 

Let us start with the trivial case of S{a,b,c,d) = 1. Linear combinations of the the selected nine which 
give the other seven vertices, have always the property that the sum of the expansion coefficients is 1. E.g., 
■~ S(a' 1) ^'-^ ~ (111: 1)7 therefore the sum of the coefficients is —'^(^a',b')ji{i,i)^'^' ~ ^^'^ 
therefore (1, 1, 1) ® (1, c', d') has the required property. Similarly 

^ a'b'c'd'{l,a',b')^{l,c',d') = {1,1,1) ^{1,1,1), (25) 

(a',6')#(l,l),(c',d')7^(l.l) 
^^'^ Z](a',b')#(l,l).(c',d')#(14) <^'b'c'd' = 1. 

Now we shall relate the above specific case with the general one. Since trivially 

E(a',6')#(i,i),(c',d')/(i,i) ^^'b'c'd'Sia', b', c', d')[S{a', b', c', a', b') ® (1, c', d')\ 

= S{1, 1, 1, 1)5(1, 1, 1, 1)(1, 1, 1) ® (1, 1, 1), (26) 

therefore the sum of the expansion coefficients for the vertex 5(1, 1, 1, 1)(1, 1, 1) ® (1,1,1) is given by 
E(a',6')5^(i,i),(c',<i')5^(i.i)"'^''^''^''^( s.m.ceY.^^^abS{a,b,c,d) = a.-n.AY,c,d'^dS{a,b,c,d) = 0, 

one has 

a'b'S{a',b',c,d) = S{l,l,c,d), (27) 

(a',6')7^(l,l) 



and 



Therefore 



^ c'd'S{a,b,c',d') = S{a,b,l,l). (28) 

(C',d')7^(l,l) 



E(a',b')#(l,l),(c',d')/(1,1) b', d, ^)S{1, 1, 1, 1) 

= -5(1, 1, 1, 1) ^(1' 1' c'' ^')c'rf' = \-S(l, 1, 1, 1)]2 = 1. 



(29) 



Similarly, one has 



= -5(a', 5', 1, l)[5(a', 6', 1, 1)(1, a', 5') ® (1, 1, 1)]. 

The sum of the expansion coefficients for the vertex 5(o', b' , 1, 1)(1, a' , b') (g) (1, 1, 1) is given by 

- E(c',d')^(i,i) c'd'5(a', b', c', d')S{a', b' , 1, 1) = 5(a', b' , 1, if = 1, 



(30) 



(31) 



i.e. again we have the required property. 

In the next step one has to show that if a{a, b, c, d) is not of the form 5(a, b, c, d) then the set of sixteen vertices 
does not belong to one face. For the assumed cr one must have either f^ab(T{a,b,c,d) ^ for some c,d, or 
Ec d c-d(j{a, b, c, d) 7^ for some a, b. Assume that the second case holds, for some a! ,b' . Then 

^ c'dV(a', 6', c',rf') 7^ CT(a', 6', 1,1) (32) 



and thus 



^ c'dV(a', b', c', d') = xa{a', b', 1, 1), (33) 

ic',d'mi,i) 
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with a; = 1 or ±3. Therefore, since one has 

E(c'.d')/(i,i) c'rf'^K: b', C, d')[a{a\ h', c', a', h') ® (1, c', d')] 
^ -a{a', y, 1, l)[cr(a', 5', 1, 1)(1, a', 6') (g) (1, 1, 1)], 

(34) 

the sum of the expansion coefficients for the vertex a{a' , b', 1, 1)(1, a', b') ® (1, 1, 1) is given by 
- E(c',d')#(i4) c'd'aia', b\ c', d')a{a', b' , 1, 1) = -xa{a' , 6', 1, l)^ = -x ^ 1, 

(35) 

i.e. the vertex does not belong to the plane defined by the chosen subset of nine! 

VI. BELL INEQUALITIES: EQUIVALENT FORM 

Above it was shown what is the full set of Bell inequalities for the considered problem. However, one must admit 
that the the derived form of the inequalities is not good for practical applications. Below it will be shown that one 
can write down the inequalities in a much simpler form. 

The full set of Bell inequalities, forming the necessary condition for the positivity of the unphysical distribution, can 
be written as 

2'' - Y.a,b.c.d=±i S{a,b,c, d)[EoQ + Eqic + Eo2d + Eiqu 
+Eiiac + Ei2ad + E2ob + E^ibc + £^22^) > 0. 

(36) 

This can also be written down as: 

2' - Ea,b.c,d=±i S[a, 6, c, d)vabcd ■E>0 (37) 

where we use the following convention 

Vabcd = (1, a, b) (g) (1, c, d) = (1, c, d, a, ac, ad, b, be, bd), 

E = {Eqo, Eqi, E02, Eio, Ell, E12, E20, E21, E22), 

and • denotes scalar product. 

We shall show that for sixteen vertices of the form S{a, b, c, d){l, a, b) (X) (1, c, d), i.e., ones defining a face, one has 
saturation of the inequality. This proves the equivalence of the inequality (j36p with the ones discussed in the previous 
chapters. Take a specific vertex belonging to the face defined by the sign function S{a, b,c,d) = X + Aa + Bb + Cc + 
Dd + Eac + Fad + Gbc + Hbd, namely 

Vaobo^o^o = S{a°, b", c", d°)il, a", b") (1, c", d") 

One has 

S{a,b,C,d)Vabcd- Vaob-c-d" 

a,fc,c,(i— ±1 

= Sia",b°,c",d") {X + Aa + Bb + Cc + Dd + Eac + Fad + Gbc + Hbd) 

a,fc,c,(i— ±1 

X (1 + a°a + b°b + c°c + d°d + a°c°ac + a°d°ad + b°c%c + b°d%d) = 2'^S{a°, b°,c°, d°) 
X (X + Aa° + Bb° + Cc° + Dd° + Ea°c° + Fa°d° + G6°c° + Hb"d°) 

= 2'iS'(a°,6°,c°,d°)2 = 2^. (38) 

That is, the inequality is saturated for all vertices belonging to the face. Since for vertices that do not belong to the 
face, that is for — Wa°h°c°!i°, one gets —2*, therefore the inequality is a tight Bell inequality associated with the face 
generated by the sign function S{a, b, c, d). QED. 
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VII. RAMIFICATIONS 

The procedure can be generahzed to an arbitrary Bell problem involving N parties, two valued observables, and 
3 local settings per observer. Generalization to more than two parties is trivial, and will be presented below. So is 
generalization to more than three settings of the derivation of the inequalities as necessary conditions for local realism 
(proof of their sufficiency may be not easy anymore). 

A. Three or more observers 

For > 2 it is straightforward to write down the full set of tight Bell inequalities in the form of a single "synthetic" 
formula. E.g., for three parties, dichotomic observables, and three settings at each side such a set of tight Bell 
inequalities is given by 

-1< J2 '?''l(a,&,c,d,e,/,5[3]) < 1, (39) 

a,b,c,d,e,f —:Ll 

where 

is any admissible "sign" function of the variables a,b, c, d, e, f with the property that in its Fourier expansion 
products ab, cd and ef never appear (i.e., all products of indices of a local observer are missing), 

• the following overcomplete bases are selected from the set of deterministic "orders" sent by the local realistic 
source to the observers: 

^[3] (a, 6, c, d, e, /)(1, a, b) (1, c, d) ® (1, e, /) = Le/,s[3i ' (40) 
where most importantly the sign function belongs to the class of admissible ones, 

• and finally 

Q[3](a, 6, c, d, e, /, ^[3]) = E ■ l^J^',^^ 5,3, , (41) 

where the set of three party correlation functions E'ooo; -E'ooi: ■■■E222 is written down as a vector (or rather a 
three index tensor, E) in R^''' = (g) R^. 

Note that the set of inequalities (f4T|) has never appeared in the literature, just like zillions of other ones derivable by 
this method. 

B. Other generaUzations 

If one wants to build Bell inequalities involving the full set of observable data for the experiment, that is in the two 
parties case, also local averages of the local results, not only the correlation functions of products results of the two 
parties, one can re-interpret the inequalities in the following way. 

Let us consider the two-observer, two-settings problem, with dichotomic observables. Consider the following subset 
of the vertices of the 3x3 problem 

W^'(m^^^TOf^;m^^^m^^^) = (1, m^^\ m^^^) ® (1, m^^\ r7?4^^). (42) 

Simply the hidden results of the local experiments for the first settings, with j = 1, 2, was replaced by 1. Now the 
averages of components of the p2)) of such vercites (or if you like convex combinations of all of them) , that is 

^P{W')W', (43) 
w 

give: normalization condition, averages of local results, and finally all correlation functions, for all combinations of the 
two pairs of settings. I.e. nothing is missing from the full description of the phenomena, for the problem in question. 
Thus we get the full description of observable phenomena for an experiment for which a local realistic model exist. 
Since the new vertices satisfy the inequalities derived for the 3x3 problem involving correlation functions, we get 
immediately a necessary condition for the existence of (I43p . In other words, we get the CH inequalities. When this 
approach is (straightforwardly) generalized to more than two parties new families of tight Bell inequalities can be 
derived. 
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C. Relation with the complete set of two settings inequalities 

The complete set of two setting Bell inequalities for N parties, dichotomic observables and two possible settings 
for each party, can be recovered from the family of inequalities presented here. Consider first the two party problem. 
One can choose the subset of admissible sign functions of the following form S{a,b,c,d) = S{a,c), i.e. depending 
on only one index per observer. Due to this fact all components of S in the expansion (|24|) which contain b and 
d drop out, and the inequality reduces to a two-setting one. Similarly, for three parties problem one can choose 
S't^l (a, 6, c, d, e, /) = S'^'^l (a, c, e), with equivalent result, etc. 
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